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PACS 42.50.Ct - Quantum description of interaction of light and matter; related experiments 
PACS 42 . 25 . Fx - Diffraction and scattering 

PACS 42 . 25 . Gy - Edge and boundary effects; reflection and refraction 

Abstract. - Despite the quantum nature of the process, collective scattering by dense cold 
samples of two-level atoms can be interpreted classically describing the sample as a macroscopic 
object with a complex refractive index. We demonstrate that resonances in Mie theory can be 
easily observable in the cooperative scattering by tuning the frequency of the incident laser field or 
the atomic number. The solution of the scattering problem is obtained for spherical atomic clouds 
who have the parabolic density characteristic of BECs, and the cooperative radiation pressure 
force calculated exhibits resonances in the cloud displacement for dense clouds. At odds from 
uniform clouds which show a complex structure including narrow peaks, these densities show 
resonances, yet only under the form of quite regular and contrasted oscillations. 



Introduction. — Mie theory is the well-known solu- 
C*-)tion of Maxwell's equations for the scattering of electro- 
magnetic radiation by spherical objects [TJ[2] . Via calcu- 
C*"Mation of the electric and magnetic fields inside and out- 
side the object the theory predicts the total optical cross 
I section, which determines the amount of scattered light, 
•and the form factor, which characterizes the far-field ra- 
diation pattern [3l[I]. Simple solutions are available in 
l regimes where the object size R differs very much from 
y—i the radiation wavelength A, or when the refractive index 
^ rn is close to unity. For example, for small phase-shifts 
. _ \m — l|i?/A <C 1 in optically dilute media |m — 1| <C 1, one 
enters the Rayleigh-Dcbyc-Gans regime whereas for small 
^ particles and small phase-shifts, one obtains Rayleigh scat- 
tering by point-like objects. 

For objects whose size is of the order of the radiation 
wavelength (e.g. water droplets in the atmosphere or in 
emulsions), Mie's full theory has to be used to find the 
scattering pattern. Mie scattering differs from Rayleigh 
scattering in several respects. While the intensity of 
Rayleigh-scattered radiation scales with the object size as 
R 6 and is identical in forward and backward direction, the 
intensity of Mie-scattered radiation is roughly independent 
of wavelength, and it is larger in forward than in backward 
direction. The greater the particle size, the more light is 
Mie-scattered into forward direction. The hallmark of Mie 



scattering, however, are the Mie resonances: those are sets 
of parameters (size, refraction index, wavelength), where 
Mie scattering is particularly strong or particularly weak. 
The sharpness of some Mie resonances make them useful 
for measuring unknown parameters such as particles' size. 

Recently, a series of papers demonstrated how collec- 
tivities of point-like Rayleigh-scattering particles may co- 
operate [5j|10| in scattering radiation into the forward di- 
rection and the relationship to Mie scattering was pointed 
out [TTJ[T2]. Here, we show that the theory of collective 
scattering by smooth distributions of point-like scattercrs 
is equivalent to Mie scattering by demonstrating that the 
premisses of both models are identical. Hence, we may ap- 
ply the Mie scattering technique to atomic clouds as long- 
as their granularity, as well as collisions and nonresonant 
atomic interactions, can be neglected. 

fn Mie theory, boundary conditions of the scattering 
object assume a fundamental role: these are generally 
sharp since dielectric spheres typically have homogeneous 
densities, while atomic clouds, in general, have parabolic 
or Gaussian density distributions and smooth boundaries. 
Within the framework of the Mie theory, we compare cal- 
culations for homogeneous and parabolic densities, and 
identify the impact of sharp or smooth boundaries on the 
occurrence and shape of Mie resonances. Our study re- 
veal that some resonances will persist, although not the 



p-1 



R. Bachclard et al. 



sharpest: this means that a jump in the refractive index 
is not required, as one could expect from classical cavi- 
ties where smooth gradients of index also allow for mode 
propagation and selection (e.g. graded- index fibers). 

Finally, we predict experimental regimes where reso- 
nances in Mic scattering can be observed in atomic clouds 
by monitoring the radiation pressure force acting on the 
cloud's center-of-mass. To this aim, we calculate the res- 
onances as a function of atom number and pump laser 
detuning. A major advantage of using resonant atoms as 
scatters is that the refraction index can be varied over 
very large ranges by changing the cloud's density and vol- 
ume, or simply by tuning the light frequency. This greatly 
facilitates the detection of the resonances. 

In this Letter, we first show the equivalence of the co- 
operative scattering by two-level atoms and the Mie scat- 
tering, assuming a scalar photon model and reducing the 
scattering equation for excitation probability amplitude to 
a differential Hclmholtz equation with complex refraction 
index. Then, we solve the Mie problem and calculate the 
radiation pressure force, before discussing the presence of 
Mie resonances in spherical clouds with parabolic densi- 
ties. Finally, we draw our conclusions. 

From collective atomic scattering to Mie scat- 
tering. — In his seminal paper [1] , Gustav Mie proposed 
an analytical solution to the scattering of light over ex- 
tended objects under the form of infinite series. This 
theory was originally developed for homogeneous media 
of refractive index mo, in which the elementary solu- 
tions to the wave equation are known to be of the form 
iprU = jn(mokor)Y n i(d, ip) (where j n is the spherical Bessel 
function, and Y n i the spherical harmonics), while outside 
they simply read j n (kor)Y n i(6,(p). The continuity of the 
fields at the boundaries eventually determines the scatter- 
ing coefficient for each (n, I) mode. 

On the other hand, the cooperative scattering by a sam- 
ple of N two- level atoms (with random position rj , transi- 
tion frequency co a and linewidth T = d 2 u!^/2TrheQC 3 , where 
d is the electric dipole matrix element), illuminated by a 
resonant uniform field is described by the following cou- 
pled equations [THS] 



d/3j 
~~dt 



r\ „ i d „ , . 



(1) 



exp(ifco|rj 



where j = 1, . . . , N and /3j is the probability amplitude of 
excitation of the jth atom, Ei(r) = _Boe lk °' r is the electric 
field of the incident laser and Ao = ujq — uj a its detuning 
with respect to the atomic transition, where ujq = cfco- In 
this approach, short-range dipole terms and polarization 
effects are neglected [13]. Neglecting granularity effects, 
the cloud can be described by a continuous field /3(r,t), 
whose steady-state regime is given by 

e lk °- r = (26 + i)j3(r) (2) 



+ | drV(r , )5B g^ (r , ) , 

where p(r) is the atomic density, 6 = Aq/T and we have 
set 

/3(r) = ^/3(r). (3) 

Let us remark that the kernel of Eq. <j2j) is the Green func- 
tion for the Hclmholtz equation, that is 



fX7 2 , ^ cxp^fcolr-r'l) 



-4ir5(r - r') (4) 



and that (V 2 + fc 2 ) exp(ik ■ r) = 0. Then, applying (V 2 + 
feg) on Eq.®, we obtain that /?(r) satisfies the Helmholtz 
equation [T4"Ill"5] 

[V 2 + fc 2 m 2 (r)]/3(r)=0, (5) 
where m(r), the cloud refractive index, is given by 

47rp(r) 



m 2 (r) = 1 



fc§ {28 + i) ' 



(6) 



Hence, the cloud of cold atoms acts on the light as a "clas- 
sical" medium of index m(r), whose imaginary part origi- 
nates in the single-atom decay term |16j : it is here respon- 
sible for the absorbing nature of the cloud, and vanishes 
only in the limit of far-detuned incident laser. 

For a cloud with spherical symmetry m(r), the so- 
lutions of the wave equation can be decomposed along 
the orthogonal basis of the spherical harmonics [1] as 
E^Lo J2's=-n u n{r)Y ns (8, tp) where the radial modes u n 
satisfy 



<(r) + 2 



u 'n( r ) 



m (r) 



n{n + I) 



0. 



(7) 



We focus on the axi-symmetric problem, where only the 
s = modes are relevant, and the solution write 



(r) = V2^TTi n (3 n u n (r)Y n0 (f). 



(8) 



n=0 



where f is a unit vector in the direction of r. As for the 
incident wave, it decomposes as 

oo 

exp^ko • r) = V^Yl V2n+li n j n (k r)Y n0 {f). (9) 

n=0 

Since the exponential kernel is diagonal in the basis of the 
spherical harmonics 



n— s— — n 



j n (k r')hh (kor) for r > r' 
j n (fc r)/i^ 1) (fc r') for r < r' 



(10) 
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where hn\z) are the spherical Hankel functions. We then 
introduce f n (r) such that 



drV(r') eX ^ |r *', l) u n (r')Y n0 (r') = f n (r)Y n0 (f) 
k a \r - r'\ 



so that, assuming p(r) = for r > R 

f n (R) = Anih^\k R) / p(r')r' 2 u n (r')j n (k r')dr'. 



(11) 



Then, the projection of Eq.© on mode n leads to 

jn(k r) = (25 + i)/3 n u n (r) + /3 n f n (r). 
Using Eq. dT^I) and defining 



A,, 



47T / p(r)r u n (r)j n (k r)dr 
>o 



we obtain 



j n (k R) 



(25 + i)u n (R)+i\ n h i n ) (k R) 



(12) 
(13) 

(14) 

(15) 



Thus, the exact solution to the scattering problem is ob- 
tained under the form of infinite series over the spherical 
modes, as is reminiscent of Mie's solution [211]. In the lat- 
ter approach, the coefficients of scattering for each mode 
are determined by using the continuity of the tangential 
electric and magnetic fields, and of their derivatives, at 
the cloud boundaries. In our problem, we considered a 
scalar electric field only, and assuming it is orthoradial, it 
reads [3] 



E s (r) = - 



^/drVM^f^^). (16) 
47re J r_r 



Now, the electric field expresses in two different ways in- 
side (r < R) and outside (r > R) the cloud 



£(° u *)(r) 



2i 
E 



^2(2n+l)i n p n X n h^(k r)P n (coim 
,. ^(2?i + l)z' l /3 n P n (cos0)4^ 



n=0 




dr'r ,2 u n (r')j n (k Q r')j h^(k r) 
dr'r'\ n (r')h^(k r')) j n (k r(jL&) 



A straightforward calculation shows that E s (r) and 
dE s (r)/dr as given by eqs. ([T7)) and (TT51) arc continuous 
at the cloud boundary r — R. This allows to conclude 
that our solution is the same as the one proposed by Mic, 
although we did not make a direct use of continuity hy- 
potheses for the electromagnetic fields. Let us also remark 
that our solution holds for any radial solution u n (r) of the 
Helmholtz equation, i.e. it applies to any spherical cloud 
with finite boundary at r = R. 



Mie resonances in nonuniform resonant media. 

— The treatment of scattering in nonuniform media is 
difficult due to the lack of explicit solutions to the wave 
equation in that case, so numerical approaches are the 
most common to tackle with this problem [3]. We shall 
here treat the case of homogeneous samples, and of those 
with a quadratic dependence of the atomic density over 
the radius. 

In the first case, the atomic density of a cloud of ra- 
dius R is po = 3N/4ttR 3 , whereas its refractive index 
is constant mo = \J\ — 3N/[(koR) 3 (25 + i)] and the so- 
lution of the wave equation inside the cloud is simply 
u n (r) = j n (komor). Thanks to properties of the Bessel 
functions Q7], A„ can be calculated explicitely as 

A„ = (2(5 + i)(k R) 2 {mojn-i(kom R)jn(koR) 

~j n -i{k R)j n (k m R)} (19) 

As for clouds of size R, with a parabolic density p(r) = 
(57V/2F)[1 - r 2 /R 2 } and volume V = 4irR 3 /3, their 
spatially-dependent refraction index reads 



m 2 q {r) 



2 i 2 2 

m c + 7 r , 



(20) 



with m c = yj\ - (15/2)A7[(fc i?) 3 (2<5 + i)\ the index in 
the core of the sample, and 7 2 = (l5/2)N/[(k R) 5 (2S+i)]. 
Remark that when k$R — > oo, but at constant central den- 
sity N/V and at fixed r, Eq.© tends toward the equation 
for the spherical Bessel function j n {kom c r): the homoge- 
neous medium limit is then recovered, as expected since 
the center of the cloud is then locally homogeneous. Using 



the substitution 



u n {r) 



- r -3/2 



w n (x), with x = /2, 



we obtain a Coulomb wave equation for w n (x), well known 
in nuclear physics [20j 



1 



2jx 



n(n + 1) 



_3_ 

To 



1 



w n (x) 



= 0. 
(21) 

Its solutions are the so-called Coulomb wave functions [5T] 
Fn/2-i/i(— m c/^l^ x )- The irregular Coulomb wave func- 
tion G is discarded because of its 1/r divergence at the 
origin, just as the spherical Hankel function is discarded 
to describe the field inside homogeneous media. Thus, we 
get the following solution for the radial component of the 
field 



47 



Let us introduce the structure factor s(k) 
direction k = fco(sin cos tp, sin sin ip, cos 9 



S (k) 



1 

TV 



dr'p(r')/3(r')e 



= s(9, ip) for a 
defined as 

(23) 



After integration over space, it turns into 

00 

s(k) = J2( 2n + l)/3„A„P„(cosfl), (24) 



n=0 
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where we have made use of 

,-2jv 







dip' sin¥ / e -*(ac OS1 /+& W) = 2irJ (yja 2 + fe 2 ) 



dy J (r'smd^l~y2)P n (y)e- i y r cos6 dy 
= 2i- n j n (r')P n (cos9). 



(25) 



Then, in the far-field limit (r> R), the scattered intensity 
simply reads 

I scat (r,6,ip) = I ^\s(e,i P )\ 2 . (26) 

Using the analytic expressions obtained for the homoge- 
neous and the parabolic density, we now turn to studying 
the resonances in Mie scattering for atomic clouds. Sim- 
ulations of the intensity scattered along the illuminating 
axis reveal that homogeneous samples exhibit some spiky 
irregular structures, whereas inhomogeneous (parabolic) 
samples are characterized by more regular oscillations (see 
FigO]). This suggests not only that the resonances do not 
only originate in the sharp boundaries of homogeneous 
obstacles, but on the contrary that they may even be eas- 
ier to observe in media with smooth densities. Actually, 
the numerous spikes in homogeneous media can be inter- 
preted as "whispering gallery modes" [22H24] , that propa- 
gate along the surface of the cloud, within a local minimum 
of the effective radial potential created by the jump of 
the medium's index. On the other hand, assimilating the 
cloud to a cavity, one can interpret the regular oscillations 
as cavity modes which survive despite the smoothening of 
the boundaries of the potential. Thus, the smoothening 
of the index leads to a modification of the local minimum 
for the effective potential which supports the whispering 
gallery modes [53], so that the latter may be attenuated 
or broadened (see FigO}, or even disappear (see Figj2]). 

Finally, the simulations confirm that the theory for con- 
tinuous media is in good agreements with TV-body samples 
(see FigUJ , although discrepancies appear for smaller val- 
ues of detuning, and thus larger refractive index. How- 
ever, simulations where N is tuned reveal that this effect 
is due to the small number of particles that one can con- 
sider numerically: since finding the stationary solution of 
(fTJ) requires the inversion of a N x N matrix problem, 
simulations are typically limited to systems of at most 
N ~ 2.10 3 particles. Thus, the finite- N effects observed 
in Figfl] are expected to vanish for larger number of par- 
ticles that one can typically find in a cold atomic cloud 
(that is N = I0 5 - 10 7 ). More specifically, one can show 
that N/a 2 ^ 1 is required to neglect granularity [511X0]. 
This highlights the relevance of the continuous approach to 
model the scattering process of large- N clouds, on which 
we shall now focus. 

Resonances in Mie scattering for cold atomic 
clouds. — One of the effects of the light scattering by 




- Quadratic density 

- Homogeneous density 
N-body (homogeneous) 



Fig. 1: Intensity scattered by the sample along the illumi- 
nating axis for homogeneous (thin black line) and quadratic 
(thick blue line) continuous densities, and for a iV-atom homo- 
geneous cloud with random atomic distribution (black dots). 
Simulations realized for N = 2.10 3 and oy = ko-sV {r 2 } = 1 
(and thus ay = ^/3/5fco-R for an homogeneous density, and 
oy = -v/3/7fco-R for a quadratic one). 



a cloud of cold atoms is the displacement of its center- 
of-mass (e.g. measured by time-of-flight imaging tech- 
niques 10, 12 ) due to the radiation pressure force, which 
can be interpreted microscopically by a sequence of ab- 
sorption and spontaneous emission cycles. More specifi- 
cally, the atom absorbs a laser photon with k = koz and re- 
emits it in a direction k = fco(sin 9 cos tp, sin9 sin ip, cos 9). 
The resulting radiation pressure force reads [5] 



tt 2 

F z = -Kk -£lm(s(k )) - 



hk n 2 N 

87rr 



f 

■ / dkcos6»|s(k) 



(27) 

where firj = dEo/h is the Rabi frequency of the incident 
field. In the framework of Mie scattering, the absorption 
part of the force reads 

= -hk ^ 5>n + l)Im(A„/U (28) 

As for the emission part, using the orthogonality of Leg- 
endre polynomials, as well as the relation 

xP n {x)P n > (x)dx = 2 — — — — , (29) 



(2n+l)(2n' + l) 



it reads 



2n 



F * e) = - hk °TN + 1 ) Rc ( A « A n+i/3«/3; + i)- (30) 



These analytical expressions then allow to investigate 
numerically the presence of resonances in Mie theory 
in large- N samples [JD1[T5], such as those produced in 
Tubingen [T^], where clouds of size ay = koy/ (r 2 ) w 10 
and N ~ 10 6 were studied. The simulations reveal that 
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the radiation pressure force actually exhibits resonances 
when the number of atoms TV is tuned, at fixed detuning 
6 (see Fig(2k). Their contrast is significantly higher at 
(negative) larger values of S yet they require larger num- 
ber of atoms. Plotting the force oscillations as a function 
of the phase shift of on-axis light rays 



(j> = k / [Re(m(0, 0, z)) - 1] 



(31) 



makes clear that <j> is the relevant variable to characterize 
the resonances in clouds with quadratic densities, i.e. cav- 
ity modes. Note that since the absorption and emission 
forces compensate at small 0, Mic oscillations actually ap- 
pear in the radiative force at values of <j> higher than for 
the scattered light (typically Air). The period of these os- 
cillations can be qualitatively retrieved using the formulae 
for the homogeneous media (|19[) . If one considers a large 
cloud k R 3> 1, then j n (x) « sin(x — mr/2)/x (for x 3> 1) 
leads to the approximate expression 



At; 



25 + i r mp + 1 

— — sm((m - l)k R) 



m 

sin((mo + l)fcoi? — nir) 



(32) 



, n 2 < k R. 



This explains why the radiation pressure force will oscillate 
with the value of the phase-shift (<fi = 2Re(mo — l)fcoi? for 
homogeneous clouds, and (4/3)Re(mo — l)koR for samples 
with quadratic densities) . Note that in the case of a dilute 
cloud such that |mo— 1| <C 1 (yet |mo—l|fcoi? can be large), 
one has <j) w -(2/3)6 c>/(4<5 2 + 1), which highlights the role 
of the optical thickness [7] 



&o 



■p- J p(0,0,z)dz 



(33) 



instead of the atomic density in the collective scatter- 
ing problem |10j . Since optical resonances will emerge 
for N/6a 2 ~ 1, the conditions for the observation reads 
TV > 5a 2 . Finally, remark that when the spatial density 
TV/of becomes significant, other effects such as near field 
dipolc dipolc interactions have to be accounted for [TTUT5] . 

Here, due to the special role of the wavelength of the 
illuminating laser on resonant atoms, the detuning 5 can 
also be used as a control parameter to generate Mie oscilla- 
tions: at odds from dielectric droplets experiments where 
it changes the effective size of the cloud k\ asei R, it here 
modifies only the index of the cloud. FigJ^i depicts this 
phenomenon when the incident wavelength is tuned, and 
one recovers Mie oscillations in that case as well. As dis- 
cussed previously, clouds with smooth (quadratic) densi- 
ties exhibit more regular and more contrasted oscillations. 
Remark that using 5 as a control parameter, the curves 
F z (<f>) or I SC at[0 = O](0) are extremely similar to those 
where TV is tuned (see Fig|2b), since only the refractive 
index is important, and the phase-shift <j) captures this 
information. 




Fig. 2: (a) Radiation pressure force as a function of the number 
of atoms TV in the sample and for different detuning, namely 
8 = —667 (black lines) and S — —50 (blue lines), (b) Radiation 
pressure force as a function of the phase-shift (f> for S — —667 
(blue curves) and square of the structure factor \s(8 — 0)| 2 
(black curves), proportional to the far- field scattered light in 
the forward direction. For both figures, the thin curves cor- 
respond to homogeneous samples, the thick ones to quadratic 
densities. Simulations realized for a r — 10. 



Finally, we conclude by pointing out that the reso- 
nances in Mie theory can be observed as well at posi- 
tive detunings: As can be observed in Figj3}D, the scat- 
tered intensity is rather symmetric in 8, which can for 
example be understood by the fact that the phase-shift 
cj) w — (2/3)boS/(4S 2 + 1) is antisymmetric in 5 (yet the 
imaginary part of the refractive index is not, thus the dif- 
ference in amplitudes). The symmetry decreases with 5, 
when the refractive index deviates significantly from unity, 
yet model |T]) then loses its validity, and a more exhaustive 
model is required to study the light scattering [T3"] . 

Conclusion. — We studied resonances in Mie theory 
in resonant media with parabolic densities, a case most 
appropriate for cold atomic clouds. It was shown that the 
Mie oscillations were more regular in the case of smooth 
(parabolic) densities, where only cavity modes survive, 
and that they can be observed as well in the scattered 
intensity as in the radiation pressure force, even though 
for the lowest phase shifts the radiation pressure force is 
a less sensitive measurement. 
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aflmic clouds. 

inally, we note that tuning the illuminating laser to- 
w rd the blue (5 > 0) can also allow for the observa- 
ti of the resonances. However, when 4irp/(25kQ) > 1, 
sample becomes strongly dissipative since Re(m 2 ) be- 
nes negative (see Eq.©): radiation trapping will then 
ei erge |25j . that will damp the scattering process and 
c( ipse the resonances. The occurence of this phenomenon 
ir cold atomic clouds will be the object of future works. 
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Fig. 3: (a) Radiation pressure force for a ay = 10 cloud as a 
function of the detuning 5 and for different numbers of parti- 
cles, namely N = 10 6 (black line) and N = 10 5 (blue line). 
The thin curves correspond to homogeneous samples, the thick 
ones to quadratic densities, (b) Intensity scattered along the 
illuminating axis \s(8 = 0)| 2 (black thin curve) as a function 
of the detuning, as well as real (plain blue line) and imagi- 
nary (dash-dotted blue line) parts of refractive index mo of the 
cloud, for a cloud of size a r = 100 and homogeneously filled 
with N = 2.10 6 particles. 



The parameters of FigsJ^t and [3^ correspond to the ex- 
perimental case of [12] , and it suggests that the resonances 
can be detected in the available range of parameters. Note 
that the main difference with this experiment is that the 
cloud was then cigar-shaped, and the density assumed to 
be Gaussian; yet parabolic densities are definitely more 
realistic than homogeneous ones to describe the Gaussian 
case. 

Furthermore, the contrast of the resonances is connected 
to the imaginary term of the cloud's index (O, as discussed 
in [15]. More specifically, the smaller the single-atom de- 
cay term is (compared to the detuning), the larger the 
contrast. Hence, large-iV and far-detuned configurations 
will generate sharper resonances, for they minimize the 
single-atom decay contribution. 

A more thorough treatment would of course require con- 
sidering vectorial fields [TS]. Nevertheless, simulations re- 
veal that the scalar approach provides an accurate descrip- 
tion of resonances in Mie scattering for large and dilute 
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